First Order Transition in the Spin Dynamics of Geometrically Frustrated Yb2Ti207. 
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Using neutron diffraction, ^^"Yb Mossbauer and /iSR spectroscopies, we have examined the py- 
rochlore Yb2Ti207 where the Yb^+ S' = 1/2 ground state has planar anisotropy. Below ~ 0.24 
K, the temperature of the known specific heat A-transition, there is no long range magnetic order. 
We show that the transition corresponds to a first order change in the fluctuation rate of the Yb^^ 
spins. Above the transition temperature, the rate, in the GHz range, follows a thermal excitation 
law, whereas below, the rate, in the MHz range, is temperature independent indicative of a quantum 
fluctuation regime. 



PACS numbers: 75.40.-s, 75.25.+Z, 76.80.+y, 76.75.+i 



Geometrically derived magnetic frustration arises 
when the spatial arrangement of the spins is such that 
it prevents the simultaneous minimisation of all inter- 
action energies . In the crystallographically ordered 
pyrochlore structure compounds R2Ti207, the rare earth 
ions (R) form a sub-lattice of corner sharing tetrahedra 
and a number of these compounds have been observed 
to exhibit frustration related behaviour |p[-p^. The low 
temperature magnetic behaviour associated with a par- 
ticular rare earth depends on the properties of the crystal 
field ground state and on the origin (exchange/dipole), 
size and sign of the inter-ionic interactions. For example, 
the recently identified spin- ice configuration [^jTojl^ has 
been linked with an Ising-like anisotropy and a net fer- 
romagnetic interaction. Most of the published work on 
the pyrochlores has concerned rare earth ions with Ising- 
like characteristics |6 -12 and there has also been some 
interest in the properties of weakly anisotropic Gd'^^ [|j . 
Less attention has been paid to the case, considered here, 
where the ion has planar anisotropy. 

To date, in systems where geometrical frustration 
may be present, two different low temperature magnetic 
ground states have been considered. First, under the in- 
fiuence of the frustration the system does not experience 
a magnetic phase transition and remains in a collective 
paramagnetic state with the spin fluctuations persisting 
as T ^ p|,^,|6| |Tl]JT^ . Second, a long range ordered 
state is reached through a phase transition which may be 
first order (|jl|,|l|l. Our resuhs for Yb2Ti207, obtained 
using neutron diffraction, ^^"Yb Mossbauer spectroscopy 
and muon spin relaxation (/iSR), evidence a novel sce- 
nario: there is a first order transition which does not 
correspond to a transition from a paramagnetic state to 
a (long or short range) magnetically ordered state. The 



transition chiefly concerns the time domain, and involves 
an abrupt slowing down of the dynamics of short range 
correlated spins; below the transition temperature, these 
spins continue to fluctuate at a temperature independent 
rate. 

We have established the background magnetic charac- 
teristics for Yb2Ti207 in a separate study |16jl^. The 
Yb'^"'" ion crystal field ground state is a very well isolated 
Kramers doublet with a planar anisotropy, g±/gz — 2.5 
where gz and g± are respectively the spectroscopic fac- 
tors along and perpendicular to a local [111] axis. The 
net inter-ionic interaction is ferromagnetic (the paramag- 
netic Curie- Weiss temperature is 0.75 K) and it is domi- 
nated by exchange (the dipole-dipole interaction is rela- 
tively small due to the modest value, 1.15/xb, see below, 
of the Yb'^^ moment). Specific heat measurements jlj] 
(Fig. |l|, top, inset) have evidenced a sharp peak (A tran- 
sition) at ^ 0.2 K with an associated magnetic entropy of 
~ 0.18i?ln2, and a broad peak centered near 2 K. The 
total magnetic entropy was estimated to be 0.97i?ln2 
|l8| suggesting there may be a small amount of missing 
low temperature magnetic entropy. 

Single phase, polycrystalline Yb2Ti207 was prepared 
by heating the constituent oxides up to 1400 °C with four 
intermediate grindings. The neutron diffraction mea- 
surements were made down to 0.11 K on the G41 line 
(wavelength: 0.2427 nm) at the Laboratoire Leon Bril- 
louin (Fig. 0) and to 0.065 K on the DIB line at the 
Institute Laue-Langevin. Fig. |^ shows that the crys- 
tal structure does not change with temperature and that 
below 0.2 K, there are no magnetic Bragg peaks, neither 
isolated nor superposed on the nuclear peaks. Had the 
moments of 1.15 fiB (to be evidenced below) undergone 
long range ordering, the peaks would have easily been 
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seen. We estimate that a magnetic Bragg peak would 
have remained visible if the correlation length had ex- 
ceeded 2 to 3 nm, so this provides an upper bound for the 
correlation length. An additional low angle contribution 
appears in the difference spectra below ~ 2 K (Fig. 
bottom). Its intensity initially grows as the temperature 
decreases but there is no significant change on crossing 
the temperature of the specific heat A-transition. If we 
interpret this additional contribution in terms of short 
range ferromagnetic correlations, then using a model of 
non-interacting spherical clusters each having a uniform 
magnetisation (solid line on Fig. 0, bottom) |19| we ob- 
tain a sphere diameter of ^ 1-5(2) nm. It is possible 
that not all of the additional signal seen for 29 < 15° 
is of magnetic origin ||2^. We consider our estimate of 
the magnetic correlation length below the specific heat 
A-transition to correspond to an upper limit. The ab- 
sence of long range order contradicts previous supposi- 
tions 
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FIG. 1. Neutron diffraction measurements for Yb2Ti207. 
Top: measured (points) and calculated (solid line) for the 
paramagnetic state at 7 K; the ticks indicate the positions of 
the nuclear Bragg peaks and the difference between the mea- 
surement and the Rietveld refinement is also shown. Bottom: 
the points correspond to the difference between the measured 
values at 0.11 K (below the specific heat A-transition) and at 
7 K. No magnetic Bragg peaks are visible. The upturn ob- 
served below 26 ~ 15° and the fitted solid line are discussed 
in the text. The inset in the top part shows the specific heat 
data taken from ref. [^. 

Selected ^™Yb Mossbauer absorption spectra are 
shown on the left panel of Fig. At 0.036 K, a five line 
spectrum is observed indicating there is a "static" hyper- 
fine field (i?hf) which we find amounts to 115 T. In the 
present case, "static" means the fluctuation frequency of 
the field is less than the lowest measurable ^"^^Yb value 
of ~ 15 MHz. Knowing that for Yb^+ the hyperfine field 
is proportional to the 4f shell magnetic moment, we ob- 



tain that each of the Yb'^"'" carries a magnetic moment of 
1.15 fiB- As there is no long range order, the hyperfine 
field is associated with the short range correlated Yb'^"'" 
moments. 

In the absence of a significant quadrupole hyperfine 
interaction, we cannot obtain the local direction of the 
Yb'^"'" magnetic moment by directly measuring the angle 
it makes with the principal axis of the electric field gradi- 
ent (a [111] direction). Instead, we make use of the prop- 
erty that for an anisotropic Kramers doublet, the size of 
the spontaneous magnetic moment is linked to the an- 
gle 9 it makes with the local symmetry axis (a [111] axis) 
through the relation: Myh = ^g±i^B/ {cos 9Vr^ + tm?9) 
where r is the anisotropy ratio gj_/gz Using the 

measured moment and the known g values, we find 
6 = 44(5)° . Thus each moment does not lie perpen- 
dicular to its local [111] axis as would be expected if 
the orientation were governed only by the crystal field 
anisotropy. The uniform tilting of each moment away 
from its local [111] axis is probably related to the com- 
bined infiuence of the crystal field anisotropy and of the 
(presumably anisotropic) exchange interaction with the 
non-coUinear neighbouring moments. 
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FIG. 2. Left panel: ^^"Yb Mossbauer absorption in 
Yb2Ti207 above, within, and below the first order transi- 
tion occurring near 0.24 K. The 7-ray energy is = 84 keV, 
and the ground and excited nuclear spin states are Ig = and 
/cx = 2, respectively (1 cm/s corresponds to 680 MHz); Right 
panel: thermal variation of the size of the Yb'^"'' magnetic 
moment obtained from the hyperfine field (top) and relative 
weight of the "static" magnetic fraction (bottom). The lines 
are eye-guides. 

When the temperature is increased to 0.23 K, an ad- 
ditional single line sub-spectrum appears. It is linked 
with the fraction of the Yb'^^ whose moments fluctuate 
"rapidly" so that the magnetic hyperfine splitting is "mo- 
tionally narrowed" . The two sub-spectra (see Fig. ^, left 
panel at 0.24 K) are both present up to 0.26 K evidenc- 
ing the coexistence of regions with "static" and "rapidly 
fluctuating" moments. Fig. ||, right panel, shows that as 
the temperature increases, there is a progressive decrease 
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in the relative weight of the "static" hyperfine field sub- 
spectrum. These features clearly evidence a first order 
transition, taking place at slightly different temperatures 
in the different parts of the sample. 

The single line (sub)spectrum progressively narrows as 
the temperature increases. As magnetic correlations are 
still present above 0.24 K, we attribute this change to 
the progressive increase in vm, the fluctuation rate of 
i/hf- The relation between the dynamic line broadening, 
AFij, and i/m, is written AF/j = {^iH^Y /vm where /i/ 
is the ^^"Yb nuclear moment As shown on Fig. 

^, when the temperature is lowered from ~ 1 K to just 
above 0.24 K, the rate decreases from ~ 15 to 2 GHz. 
This decrease is linked to the slowing down of the spin 
fluctuations which accompanies the development of the 
short range correlations. Below 0.24 K, vm drops to a 
value which is less than the lowest measurable ^^°Yb 
Mossbauer value. 

The /^SR study (23) was carried out at the ISIS facility 
(some additional measurements were also made at the 
PSI facility) over the range 300 to 0.04 K, mostly in a 
longitudinal field of 2 mT so as to decouple the small con- 
tribution of the nuclear spins to the measured depolarisa- 
tion. Fig. ^ shows typical results for the time dependence 
of the asymmetry which is written aP°''P(i) = azPz{t) + 
Obg where the first term on the right hand side originates 
from the sample and the second is a temperature inde- 
pendent constant background contribution linked with 
the muons stopping in the silver sample holder and in the 
cryostat (obg — 0.065 for ISIS and close to for PSI). 
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FIG. 3. Panel a: /xSR spectra recorded at ISIS in a lon- 
gitudinal field of 2 mT. A marked change occurs on crossing 
the temperature (~ 0.24 K) of the specific heat A transition. 
Panel b: short time part of the PSI data at 0.200 K confirm- 
ing there are no short time oscillations. The slight difi'erence 
visible between the ISIS and PSI data at 0.200 K is linked 
with the first order nature of the transition and the different 
thermal and magnetic field histories of the two experiments. 
The dashed and solid line fits are described in the text. 

From 300 to 0.275 K, Pz{t) is well represented by an 
exponential relaxation function (solid lines at 4.3 and 
0.275 K on Fig. |a): Pz{t) = exp(-A^i) where is 



the muon spin-lattice relaxation rate. From 100 to 
4 K, Xz is small (~ 0.1 MHz) and in keeping with the 
paramagnetic nature of the Yb^+ , it shows little thermal 
dependence. As the temperature is lowered below ~ 1 K, 
\z increases progressively to reach 0.52(3) MHz at 0.275 
K. The origin of this increase is the same as that for the 
dynamic broadening AF^ of the Mossbauer spectra, i.e. 
the slowing down of the electronic spin fiuctuations. 

As the temperature is lowered below 0.275 K, Pz{t) 
first becomes moderately non-exponential and then 
towards 0.2 K it abruptly becomes strongly non- 
exponential (Fig. ^). Below 0.2 K, Pz{t) is independent 
of temperature and there are no spontaneous oscillations 
in accordance with the neutron result that there is no 
long range order. The spectrum of Fig. ^b which extends 
to shorter times, clearly confirms this result. The shape 
of Pz{t) below 0.2 K is reminiscent of a dynamic Kubo- 
Toyabe (KT) decay |Q associated with a slowly fluc- 
tuating ensemble of moments, the decrease of Pz{t) be- 
yond ^ 1 /iS indicating that fluctuations are still present. 
The KT decay provides an approximate account of the 
shape of Pz{t) but there is a noticeable misfit below 0.5 
/iS (dashed lines on Fig. H). A similar misfit was also 
observed for Y2M02O7 near 0.03 ^s at 2.5 K (l3|. A bet- 
ter fit is provided by the Gaussian Broadened Gaussian 
(GBG) model |^^, where the single Gaussian distribu- 
tion of the KT model is replaced by a collection of dis- 
tributions (solid lines at 0.200 K on Fig. ||). From the 
fit of Fig. the mean value of the GBG distribution 
is Alt/t^i — 5.7 mT, the ratio of the width of the col- 
lection of distributions to the mean value is ~ 0.38 and 
the electronic fluctuation rate is ~ 1 MHz. All are 
independent of temperature below 0.2 K. 

To obtain an estimate of the effective value of the spin 
fluctuation rate, u^, at temperatures above 0.24 K, we 
insert the experimental values for (obtained in 2 mT) 
into the expression A^ = 2A^rj,/;y^. The value for 
then depends on the value of Aht, the Yb'^+ - muon 
spin coupling above 0.24 K. With the choice Aht/t^i 
= 31.8 mT, the obtained from the muon measure- 
ments scale remarkably well with the vyi obtained from 
the Mossbauer measurements (Fig. Their common 
thermal dependence is well fitted (solid line on Fig. ^) 
by the thermal excitation law v = I'o exp[— £'bh/(^Br)], 
with j/Q = 17 GHz and mean barrier height Ebh = 0.5 K. 
The associated drop in at 0.24 K amounts to a fac- 
tor of lO'^. Choosing a different value for Aht/T/^ would 
lead to different rates above 0.24 K and to a different size 
for the drop. For example, with Aht/7/j = 5.7 mT (the 
value we obtained for Alt/7/j below 0.24 K), the drop 
is ~ 10^ (this is also the minimum size of the drop al- 
lowed by the Mossbauer results) and with Aht/7/j — 85.6 
mT, obtained by scahng from TbaTiaOy ||, it is ~ lO'^. 
Combining the Mossbauer and the muon results thus in- 
dicates that the spin fiuctuation rate undergoes a first 
order change by a factor of 10^ to 10'^. We mention that 
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between 2 and 0.275 K, we find that \z does not have the 
usual Lorentzian dependence on the apphed longitudinal 
field. This anomalous behaviour will be discussed in a 
separate publication. 
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FIG. 4. Estimate of the Yb^"*" fluctuation rates as obtained 
from ^'^"Yb Mossbauer {vm) and /iSR {vfj) measurements. 
The first order change in the fiuctuation rate takes place at 
the specific heat A transition. Below ~ 0.24 K, the fiuctuation 
rate is independent of temperature and has dropped below the 
lowest value which is measurable with the Mossbauer method 
(dashed line). The solid line follows a thermal excitation law 
as described in the text. 

The low temperature magnetic properties of Yb2Ti2 07 
are therefore unusual. In the short range correlated re- 
gion from 2K to ~ 0.24K, the spin fluctuation rates 
follow a thermal excitation law. This crystallographically 
ordered compound thus possesses barriers against spin re- 
orientation. On crossing the temperature of the specific 
heat A-peak [Q, no long range order appears, but there 
is a first order drop, of two to four orders of magnitude, 
in the fluctuation rates of the correlated moments. As 
T — > 0, the fluctuations persist at a temperature inde- 
pendent rate of ~ 1 MHz, and they take place between 
directions which make an angle of ~ 44° relative to the 
local anisotropy axis (a [111] direction). The fluctua- 
tions of the moments thus conserve this angle and hence 
involve spin flips or spin spirals around a [111] axis. 

Above 0.24 K, we speculate that the spin fluctuations 
could link more random directions. The flrst order drop 
in the fluctuation rate would then be associated with a 
change in the nature of the fluctuations such that below 
0.24 K, the system explores a subset of the states that are 
explored at higher temperatures. The described transi- 
tion, first order without the appearance of any long range 
correlations and with a change of dynamics, evidences 
properties which parallel aspects of the general gas-liquid 
transition or of the specific water-ice transition. To our 
knowledge, such a transition has not previously been ob- 
served in a magnetic system. 

When exotic low temperature properties have previ- 
ously been observed in the rare earth pyrochlores, they 
have related to ions with Ising character (Tb^+, Dy^+, 



Ho'^"'"). The present results show that further novel be- 
haviour is evidenced when the rare earth has a pla- 
nar anisotropy. Although Yb2Ti207 does not corre- 
spond to the spin-ice scenario (in addition to the pla- 
nar anisotropy, there is, at most, only a small missing 
entropy), it shares with the spin- ice state the character- 
istic of short range correlated moments which continue 
to fluctuate between specific directions as T ^ 0. Here 
however, the low temperature state is reached through a 
novel route: a first order dynamical transition which can 
be viewed as a change from a thermally excited regime 
to a low temperature quantum (or tunneling) regime. 

We thank Isabelle Mirebeau and Collin Broholm for 
useful discussions and Nadine Genand-Riondet for tech- 
nical assistance. 
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